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ABSTRACT 


Many evolutionary steps are required before a solar-thermal power 
system can be moved from the research laboratory phase to the field 
applications area. It Is the objective of this paper to discuss some 
possible means for achieving the required technology. Simplifications 
In construction techniques as well as In measurement techniques for 
parabolic trough collectors are described. Actual measurement data 
Is also given. 
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I. INTRODUCTION 


The objective of this peper Is to study the performance of a parabolic 
trough solar- thermal collector with axis oriented parallel to the N-S Earth axis. 
The emphasis Is on simplified experimental techniques to provide rapid evaluation 
of a collector system. A brief description of a simple parabolic collector Is 
given In Section II. In Section III a simple method for determining the thermal 
losses In a receiver pipe Is described. 

An analytical formulation for the direct component of Insolation on a clear 
day Is given In Section IV In preparation for the actual determination of the 
performance of a solar-thermal collector, which is given In Section V. In 
Section VI actual experimental data Is analyzed by the methods developed In the 
paper. 


II. PARABOLIC COLLECTOR DESIGN 

The basic approach In the fabrication of a parabolic collector was to 
machine parabolic ribs In a programmed milling machine; these ribs were then 
attached to a torque tube as shown In Figure 1. Aluminum sheets 3.2 iin (1/8 In.) 
thick were then rolled to the appropriate circular radius to match the parabolic 
shape of the ribs. The preformed aluminum sheets were then bolted on to the ribs 
to serve as the foundation for the parabolic collector. The mirror (not shown) 
was a 3.2-mm (l/8-1n.)-th1ck sheet of acrylic aluminized on one side (Ram Products: 
Industrial grade mirror) and was clamped Into the structure In a simple manner for 
easy replacement. The aluminized acrylic mirror had a reflectivity of around 85* 
even though the radiation had to traverse the plastic twice. The protection of the 
mirror by the plastic Is a desirable feature for hostile environments and well worth 
the small sacrifice In reflectivity. 
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The drive mechanism was not a matter of concern in this investigation, and 
a synchronous motor with a gear box was used in conjunction with a chain drive 
to provide adequate solar tracking. The techniques useo profited from the 
experiences of the Minneapolis -Honeywell project^^^*, and the author is grateful 
to Dr. J. Ramsey for helpful information. 

The primary objective of this project was to develop measurement techniques, 
and the collector was designed to facilitate these measurements. No effort 
was made to optimire the fabrication procedure described. Even so, the performance 
of the collector met all expectations as will be described. 

HI. DETERMINATION OF THERMAL LOSSES 

A simple way to determine the thermal losses in a solar-thermal collector 
is shown in Figure 2 . A thin-walled stainless steel tube is treated with the 
appropriate selective coating and mounted as shown. A high-current low-voltage 
power transformer is then used to heat the tube. After temperature equilibrium 
is reached, the thermal loss is measured directly by measuring the electrical 
power being supplied. Thermocouples are placed inside the tube at a few 
places to ensure that uniform equilibrium temperature has been achieved; this 
is necessary to ensure that no longitudinal heat flow is occurring. Figure 3 
shows a typical set of curves for a 2.5-cm stainless steel tube prior to 
application of any coating. Figure 4 is a photograph of the experimental 
arrangement used. 

Lest one dismisses the measurement as being trivial, it should be pointed 
out that the measurement of the effective infra-red emissivity of a selectively 
coated surface is not an easy task. With the technique described here it is 


Numerals in parenthesis refer to references cited at end of paper. 
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only necessary to evacuate the enclosing glass tube and to know the enissivlty 

( 2 ) 

of the glass' ' for the corresponding tenperature. One can then cocrpute the 

emissivlty sought. Another difficult measuretnent which can be perform, ed quite 

easily with the apparatus described Is the behavior of convection loss as a 

function of tilt angle relative to the horizontal position. It should be 

appreciated that these measurements are difficult to perform once the collector 
pipe Is mounted on the parabolic collector assembly. 

Finally, It Is necessary to point out that there will be a small difference 
In the thermal properties of the collector pipe as measured here compared with 
the actual Installed pipe. Here the temperature Is uniform In the circumfer- 
ential d^rect1on at any given location of the pipe. In the Installed pipe, 
the temperature on the side where the sunlight Is focussed may be higher than 
on the opposite side. However, the temperature difference Is not large enough 
to yield serious errors. 

IV. ANALYTICAL REPRESENTATION OF INSOLATION 
In preparation for the measurement of the performance of a solar energy 
collector, It Is desirable to derive a simple analytical expression for the 
clear day Insolation. For concentrating collectors only the direct normal 
component of the radiation Is of Importance. 

The mathematical expression for the direct normal solar Intensity as given 
In the ASHRAE Handbook of Fundamentals' ' Is compact but cumbersome for calcu- 
lations. In addition It Is Inaccurate at sunrise and sunset. This becomes 
obvious when one plots 

Ijjl^ ■ A exp (-B/s1n b) (1) 

where 

A ■ apparent solar Irradiation at (at air mass ■ 0} 

B ■ atmospheric extinction coefficient 
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and 


s1n b ■ cos L cos6 cos H + sin L s1n6 
here L ■ local latitude 

6 ■ solar declination 

H ■ hour angle ■d/4)x no. m1n from local solar noon. 
Equation (1) Is based on a flat Earth-atmosphere approximation, and Is 
at least useful over the Important part of the day. A more accurate repre- 
sentation is easily derivable for a uniform, concentric, spherical atmosphere 
of finite thickness. This Is given by 

" A exp (-B'S) (2) 

where 

S - ♦ ^1 + ^ cot^ 0 -1 J sine (3) 

A ■ same as ASHRAE 

R * Earth's radius (64*^0 km (4000 ml)) 
h • thickness of a uniform atmosphere (32 km (20 ml)) 

0 > solar angle measured from local horizon 
and 

B'h - ASHRAE B 

Unfortunately, even If Eq. (2) Is more accurate than Eq. (1), It Is evident 
that neither are suitable for computational purposes. Fortunately, a simple 
expression can be used as a good approximation to Eq. (2). This Is given by 

’dn ■*'[’- (-8" •>] 

0 <9 < i mornings 

r 1 

IpN ' A* [l - exp (-B" (tt - 0 ))J 

J < 0<r afternoons 
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where the constants A* and b" are adjusted to match the data on Insolation. 
Figure 5 is an actual recording from a tracking heliometer and shows that the 
direct-normal component of insolation it; indeed well represented by Eq. (4). In 
Equation (4) the angle 0 is given by 


0 - (0 t 
0 


(5) 


where w is the Earth's rotational rate which is, of course, a seasonal variable. 
0 

Equation (4) will be used in calculating the performance of a solar energy 
collector. The utility of Eq.(4)may be illustrated, for example, in computing 
the total solar energy available per annum. This is given by 


W = 2 


360 


'Ll 

n =1 


- o)n 

n n 




( 6 ) 


where a' = A„ + a. cos n (7) 

where the constants A^ and a^ may be computed from ASHRAE data, and n is the day 
number in degrees with Dec. 21 taken as day zero and day 360. (The fact that 
there are really 365 days per year is a minor defect in the Earth's orbital 
motion which we overlook in the interest of computational elegance.) Similarly, 
b" is given by 

b" * - bp cos n (8) 

Again the constants B^ and b^ may be computed from ASHRAE. The calculations for 
Wn and t^ are likewise very simple. Figure 6 shows a plot of the ASHRAE constants 
A, B, ard 6, the tilt angle of the Earth. 
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The Integration in Eq. (6) Is easy; this is not the case for either expression 
in Eq. (1) or (2). 


V. SOLAR COLLECTOR EVALUATION 

The performance of solar energy collectors may be determined through some 
simple calorimetric measurements. By supplementing experimental data with 
reasonable approximations it is possible to derive a theoretical model which 
permits the calculation of the performance under various operating conditions. 

The analyses described here are particularly adapted to cylindrical trough 
collectors. As stated p-'evfously the detailed design and fabrication of solar 
collectors have been discussed in considerable detail in the literature (Refs. 1,4). 
The purpose here is to consider the evaluation of the performance of a given 
solar-thermal system. A simplified calorimetric procedure is suggested as shown 
in Figure 7. Instead of using the usual flowing fluid for measurements, a 
solid rod of metal, such as aluminum, is used for calorimetry. 

The thermal balance equation for Figure 7 is given by 

^ * ^1 ( ^ ^2 (t - T,) (5) 

where 

I = direct normal insolation 
0 

Ag = effective area of solar collector per unit length 

rig = optical efficiency 

M = mass of A1 rod 
C = specific heat of A1 rod 
0 • solar angle (see Figure 7) 
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T ■ tenp of od in Kelvins 
T, ■ ambient temperature 

O 

K-| s radiation loss parameter 

Kg = convection loss parameter (in linear approximation) 

The motivation for the present approach is the observation that the inso- 
lation is substantially constant between 10 a.m. and 2 p.m. on a clear day. 

For a polar-mounted collector, the solar angle, 0, is also constant; hence, 

the input is an easily measured constant for any given clear day. Under 

these circumstances the non-linear heat balance equation may be solved in 

various approximations as discussed below. 

A. THE LINEAR APPROXIMATION 

It is well known that for small variations in temperature the radiation 
loss term in Eq. (9) may be approximated as 

'<1 f’o) 

where T^ is the mean operating temperature. Thus for small variations in T 
about T^, the radiation loss (10) may be combined with the linear loss term 
in Eq. (9) or 

I(t) = HC ^ + K (T -Tj) (11) 

where K is a new constant which accounts for all losses and 

Although Eq. (11) is applicable for a relatively small temperature range, it is 
exactly solvable even for the case when insolation I(t) is a variable (cloudy 
day). Some insight into the dynamics of a solar collector could thus be gained. 
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B. NON-LINEAR APPROXIMATION 

An alternative approximation to Eq. (9) is to assume dominance of the 
radiation loss term; then 

I(t)-HC^ *Ko (12) 

where is adjusted to account for all thermal losses. Eq. (12) is applicable 
especially to high-temperature operation of a collector. Figure 8 shows a 
typical heating curve for constant insolation. The first deduction from 
Fiqure 8 and Eq. (12) is that initially there are no losses because T * T,. 

Q 

Thus 



= I^A^n^ cos 0 


(13) 


where (dT/dt)^ indicates the initial slope at t = 0. All quantities in Eq. (13) 
are known except the optical efficiency, which is given by 



here 

p = mean reflectivity of mirror 
T = mean transmissivity of glass tube (see Figure 7) 
a = mean absorptivity of pipe 

The second observation from Figure 8 and Eq. (12) is that when thermal 
equilibrium is achieved, all the solar pov/er incident on the receiver pipe is 
reradiated as infrared power and 
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( 15 ) 


where K 1s the hitherto unknown loss parameter, and T Is the stagnation temper- 

0 i 

ature. All quantities in Eq. (12) are now known, and for any operatinq temper- 
ature, T^, the thermal power available is given by 

U-^a) 

In these equations it is to be noted that really stands for the apparent 
ambient temperature. Unfortunately this effective temperature is not directly 
measurable. It is a function of the ambient temperature and the temperature 
of the glass tube, which varies from point to point. The infra-red properties 
of the glass finally provide the weighting function to be used with the above 
two temperatures. 

One way to obtain, experimentally, a mean value for T_ is to record the 

a 

cooling curve as well as the heating curve. This is simply achieved by turning 
the solar collector 180 degrees away from the sun after the rod has reached 
stagnation temperature. 

The interesting information in the temperature vs. time recording is, of 
course, in the high-temperature region. Equation (12) is exactly solvable for 
high temperatures, and the solutions are, for heating and cooling, respectively 


t 


T = T^ tanh 
s 

L-Rt — 

+ tanh ^ ^ , 

(17) 


and 

T = T, coth 

a 

ITW 

+ coth ^ 1. 

a 

(18) 


The time t for the cooling curve should, of course, start at zero when 
T * Tg at which time the collector is turned away from the sun. 
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Summarizing this lengthy discussion. It Is observed that the heating and 
cooling curves are used with Eqs. (15), (17), and (18) to derive values for 
K and T such that the curves are well matched by Eqs. (17) and (18) at 

0 o 

least for the higher temperatures. Then Eq. (16) may be used to compute the 
available thermal power at any temperature. An actual recording from an 
experiment will be analyzed In detail in the next section to Illustrate the 
procedure. 

In the discussion of the linear approximation It was pointed out that 
Eq. (14) was exactly solvable even when the Insolation varied with time. The 
present discussion concludes with the observation that Eq. (15) Is solvable. In 
the high- temperature approximation, by rewriting It In the approximate form 

I(t) - HC - T^) (19) 

where T 1s the mean collector temperature. Eq.(19)1s reducible to the standard 
Ricatti form and Is solvable by standard methods even when the insolation varies 
with time. 


VI. EXPERIMENTAL RESULTS 

Figure 9 shows an actual heating and cooling curve obtained on a hazy day. 
The Important parameters were as follows: 

Date: June 24, 1975 
Site: JPL, Pasadena, Ca. 

2 

I = Average Insolation: 740 watts/m 
° 10 a.m. to 2 p.m. 

2 

A^ = 2.3 m = effective area of reflector 
e 

cos (9 = .92 
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Tg - 120°C (approximate plass temperature) 

T, - 30° C 

d 

C « Sp. heat of A1 » .25 Cal/gm °C 

L ■ 2.3 m * length of A1 rod 

d ■ 3.8 cm * diam. of A1 rod 

M « 7000 gms. 

The calculations proceed very simply as follows: 

(a) The total radiation Impinging on the reflector surface Is 

cos 0 « 1566 w. 

0 e 


The net power which Is converted to heating the rod Is found by measuring the 
initial slope (AT/At)^1n the recording: 


P^-HC 


(StV 


1018 watts 


This gives the optical efficiency n^for the collector as 

. _ ’’t _ 1018 _ 

% ■ 1 A cos 7 • TO 
0 e 


Similarly, the available thermal power for any temperature T^ Is found from the 


recording by measuring the slope (AT/At) 


T »T. 


or 


Alternatively, It Is possible to use Eq. (16) by first computing the value for 
Kq by use of Eq. (15) 


JPL Technical Memorandum 33-748 


n 


1018 


^oVo ** 



5.1 X 10 


-9 


673^ - 303" 


Note that If 1s significantly higher than T , the constant K Is not 
too sensitive to errors In Table 1 shows values for as calculated by 
the two methods described. It Is noted that for greater accuracy the effective 
ambient temperature ought to be treated as a variable, but the added complexity 
would not warrant the exercise. 

Figure 10 shows the available thermal power for the two operating tempera- 
tures used In the preceding examples. The choice of 230°C as a possible working 
temperature for steam was selected as a reasonably low temperature while providing 
a high enough pressure for use with a steam turbine. 

Finally, Figure 11a shows a drift curve for the collector. This curve Is 
derived as follows: the parabolic collector Is fixed In position near solar 
noon; the sun is then allowed to drift through the focusing position to provide 
the heating curve of Figure 11a. As a rough approximation the energy Impinging 
on the ro^ Is given by 


I(t) « I^A^ngCos j5 sin u)t (20) 

where tot is the angular drift rate of the sun (15° per hour). Since the temperature 
of the collector system doesn’t get very high, the rate of heating of the rod Is 
given by: 

FT Vo ^ 

Therefore the temperature as a function of time Is 

where T„^ Is the maximum temperature in Figure 11a. From Eq. (22) it Is seen 

iTIoX 

that the derivative of the drift curve gives the sinusoidal response of Eq. (21) 
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and Is shown 1n Figure 11b. The width at the half power points Is a rough 
measure of the accuracy required from the solar tracking mechanism. Here It Is 
seen to be 1 . 5 degrees . 


VII. CONCLUSIONS 

The real significance of this Investigation Is that a solar-thermal collector 
has been used as an analog computer to solve a difficult heat transfer problem. 

The rigorous solution of the heat balance equations, which treat the effective 
ambient temperature as a variable, Is not too Important because the collector 
system analyzed Is only an approximation to the real working system. 

The emphasis here has been to treat the collector as a unit and to derive 
reasonably accurate measurement techniques. The details of mirror surface 
tolerances, receiver pipe coatings, anti reflective surfaces ”For the glass tube, 
etc., have been discussed elsewhere In great detail and are not repeated here. 

At the time when the experiment was performed to obtain Figure 9, It was 
realized that the black chrome surface could be damaged by the high temperature 
of 400°C which had been achieved. Subsequently during attempts to repeat the 
experiment to obtain comparative data, for example, for the no vacuum mode of 
operation. It was noticed that considerable deterioration had occurred to the 
black chrome. 
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Table 1. Available thermal power at two operating temperatures as calcu- 
lated from graph of Figure 9 compared with calculated values 
from Eq. (16). The effective ambient temperature Is an unknown 
quantity somewhere between 30® and 120®C. 


Available 

Calculated 
from 
Fig. 9 

Calculated from Eq. 

06) 

power 

T, - 30®C 
a 

T. - 60®C 

fl 

T. » 90®C 

A 

P^(230®C) 

733 W 

728 W 

749 W 

776 W 

P^(300®C) 

554 W 

500 W 

521 W 

548 W 
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VACUUM 



Fig. 2. Method for measuring thermal 
lo'.ises in a solar-thermal collector 



TEmRATURt Of PIPE, *C 

Fig. 3. Typical thermal loss curves for a 
stainless steel pipe without surface 
treatment (heat loss per unit length of 
2.5-cm-diam stainless steel pipe) 
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VACUUM 

OAUOf 


MC I, 1974 I 
GOIDSTONI, CALI7 


• POINTS USID TO 

CAICULATI CONSTANTS 
IN EQ. (4) 

O POINTS USED TO CHECK 
EQ. (4) AGAINST DATA 


Fig. 4. Experiment*) 1 arrangement for measuring heat losses 
In a collector pipe 


AMMITU 


HIOH<U*tfNT 
POWtH SUPPLY 


PACIPIC STANDARD TIME 


Fig. 5. Typical clear day direct normal 
Insolation at Goldstone, California 
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MONTH (2IST DAY Of f ACH MONTH) 


Fig. 6. Plot of ASHRAE constant A, B. ^ (Note that A has been 
converted to kw/m^) 



Fig. 7. Test arrangement for a 
solar- thermal collector 
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Fig. 8. Hypothetical heating curve for 
arrangement of Figure 7 



HEATING », min COOLING mm 

Fig. 9. Heating-cooling curve? for A1 rod 
coated with Harshaw Chemical "black 
chrome" 
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Ig, W/m‘ 


Fig. 10. Available thermal power as a 
function of insolation calculated from 


Eq. (16) 



Fig. 11. (a) Solar drift curve (see text); 

(b) derivative of drift curve for deter- 
mining tracking accuracy required (executed 
freehand to illustrate a procedure) 
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